MERGING RATES OF COMPACT BINARIES IN THE UNIVERSE: 

ON : GRAVITATIONAL WAVES AND GAMMA-RAY BURSTS 
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Merging rates of compact binaries (double neutron stars or black holes) are calculated based on 
the modern concept of binary stellar evolution. It is found that the initial laser interferometers 
^ with an rms-sensitivity of 10~^^ at the frequency 100 Hz can detect 10-700 black holes and only 

, ~1 neutron star coalescences in a 1-year integration time. Implications of the evolutionary 

effects to the cosmological origin of GRB are also discussed. 



H ' 1 Introduction 



In a few years several initial ground-based laser interferometers aimed at searching for gravita- 
tional waves (GW) will start working (LIGO0, VIRGol, GEO-60oi, TAMA-30oEI), so at present 
the question about what kind of events and how frequently will the interferometer register is 
very important. 

Undoubtedly, the most reliable GW sources are merging compact binary stars - double 
neutron stars (NS) and black holes (BH) of different stellar masses. In the same time, merging 
compact binaries may underly the origin of cosmic gamma-ray bursts (GRB) (Blinnikov et a/.H, 
Pazcyhski§, Meszaros and Rees@). 

In Section ^ we briefly describe observational and theoretical data about double degenerate 
compact binaries and estimates of their galactic merging rates. In Section |^ we discuss evolu- 
tionary scenario parameters that mostly affect binary NS and BH merging rates. The results 
of calculations of the merging rates in our Galaxy are presented in Section ^. Then we discuss 
the transition from the merging rates in an individual galaxy to the detection rate by a GW 
detector with characteristics similar to those of the initial LIGO/VIRGO interferometer (this 
is the most important result of the paper about GW events). Two last parts of the paper deal 
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with GRB (Section |5|) and unidentified soft X-ray flashes (Section ^) recently reported in the 
Einstein IPC archive studies!. 

2 The Present Situation in the Galaxy 

2.1 Observational Data 

At present, the following facts are known about double degenerate compact binaries in our 
Galaxy: 



Table 1: Binary PSR with NS secondaries (Nice et al.U) 
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Table 2: BH Candidates (CherepashchukB). 
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1. A few binary radiopulsars are known to have the secondary NS component (Table [T|). 

2. Three of these binary pulsars must coalesce due to the orbital angular momentum removal 
by GW on a time scale shorter than the age of the Universe (the Hubble time tn — 
15 • 10^ yr). 

3. No binary pulsars with BH is known as yet (although from evolutionary considerations one 
may expect one such an object to be formed in the Galaxy per about 1000 single pulsars, 
Lipunov et a/. 1111) 



4. 



No binary BH has been found so far. 



5. In contrast, 10 BH candidates are already known in X-ray binary systems with normal 
companions 0. Their parameters are listed in Table Note that the mean BH mass in 
these systems is <Mi,h>— 8.5Mq , i.e. BH formed in stellar evolution are notably more 
massive than NS (with the typical mass I.4M0). 

2.2 Binary NS Merging Rate Estimates 

At present, it is possible to estimate binary NS merging rate in two ways: using the binary 
radiopulsar statistics observed and making various computations of binary stellar evolution. 

"Observational" estimates 



Phinneylll . . . . 1/10*' yr 

Narayan et alM 1/10^ yr 

Curran and Lorimer0 ^ 3/10^ yr 

van den Heuvel and Lorimerll3 8/10^ yr 

"Bailes limit" 111 < 1/10^ yr 

"Theoretical" estimates 

Clark et alM 1/lOM/lO^ yr 

Lipunov a/. Ha 1/10^ yr 

Hils et alM 1/10^ yr 

Tutukov and Yungelson@ 3/10''-l/10'^ yr 

Lipunov et alM < yr 

Portegies Zwart and SpreeuwH 3/10^ yr 

Lipunov et alM 3/10^-3/10^ yr 



We emphasize that although theoretical merging rates are systematically higher than ob- 
servational ones, both estimates do not contradict each other. The main argument is that the 
first (observational) estimates of binary NS merging rate are based on the statistics of binary 
systems, in which only one of the components shines as radiopulsar, which is not at all the nec- 
essary condition for merging to occur. Thus the observational estimates are in fact only lower 
limits. There are also some selection effects that can change estimates from both groups. 

3 Parameters of binary compact star evolution 

To calculate binary evolution, we have used the population synthesis method (the Scenario 
Machine code), which is in fact a version of the Monte-Carlo calculations. Here we shall not 
enter into detail of the evolutionary scenario used. Much more detailed description of the method 
can be found in our reviewed. 

An example of the evolutionary track leading to BH-I-BH binary system formation is shown in 
Fig. ||. A short glance to this track is sufficient to understand that there are a lot of evolutionary 
scenario parameters, which affect different stages of the binary evolution. Fortunately, a very 
limited number of parameters has effect on the compact binary merging rate. 

The most important (and practically unique) parameter changing the galactic binary NS 
merging rate is the distribution of an additional (kick) velocity imparted to NS at birth. The 
kick velocity distribution widely accepted at present is derived from the analysis of spatial 
velocities of single radiopulsarscB. 

We have approximated this 3-dimensional distribution as 

fLL{x)dx OC x0-19(l + x6-^2)-l/2^^ 
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Figure 1: Massive binary evolutionary track. 



where x = w/wq and the characteristic velocity wq is a parameter in our calculations. The 
observed pulsar transverse velocity distribution corresponds to wq = 400 km/s. 

In contrast, for BH two additional parameters appear. First of them is a threshold main 
sequence stellar mass Mcr for the star to collapse into a BH after its nuclear evolution has 



ended. This parameter is sti 



poorly determined and varies in a wide range: e.g., according to 



van den Heuvel and HabetsE3 Mcr = 4O-8OM0; Tsujimoto et a/. EH give 4O~6OM0; Portegies 



Zwart et alM denve >2OM0. 

The second parameter is the fraction of the presupernova mass, ksH, collapsing into BH. 
This parameter is fully unknown, so we varied it from 0.1 to 1 in our calculations. 

The third parameter, as for NS, is the kick velocity. Clearly, in the general case the more 
massive BH will acquire smaller velocities than NS. In our calculation we used the following ad 
hoc relationship 

WBH _ MpreSN " MbH 
WnS MpreSN - Mqv ' 

where Mqv = S.SMq is the maximal NS mass (Oppenheimer-Volkoff limit). When BH mass is 
close to Moy, velocities of BH and NS are assumed to be almost the same, whereas at kBH = 1 
BH kick velocity is assumed to vanish. (Of course, other dependences wbh/wns are possible, 
but their specific shape weakly affects the results). 



4 Detection rate of binary compact star merging 



Under the assumptions made above, we can calculate the binary merging rates in the Galaxy 
R. The results are presented in Fig. ^ and ^ 




Age, yr 



Figure 2: The dependence of different compact binary merging rates in a model "elliptical" galaxy (where all the 
stars were formed simultaneously at the moment t = 0) with mass M = W^^Mq. The horizontal line shows the 
mean merging rate of binary NS in a spiral galaxy of the same mass with a stationary star formation. 




Figure 3: The dependence of different compact binary merging rates in a spiral galaxy with 10 M0 on the 

characteristic kick velocity wo- 




Figure 4: The total merging rate of NS+NS, NS+BH, and BH+BH binaries as would be detected by a laser 
interferometer with hrms = lO"^'^ for Lyne-Lorimer kick velocity distribution with wq — 200-400 km/s and 
BH progenitor's masses M, = 15-50Mq, for different scenarios of binary star evolution as a function of kth- 
NS+NS mergings are shown separately. In all cases BH+BH mergings contribute more than 80% to the total 
rate. The filled "Loch-Ness-monster-head"-like region corresponds to BH formation parameters M» > ISM© 

and = 0.5. 



After having found the merging rates i? in a typical galaxv. we need to go over the event rate 
D at the detector. Applying the optimal filtering technique ej, the signal-to-noise ratio S/N at 
the spiral-in stage is 

S M^J^ 
— oc — - — . 

N d 

Here M^h = (Mi Ms (Mi + Ma)^/^ is "chorp" -mass of the binary system. This means that 
for a given S/N our detector can register more massive BH from larger distances than NS. The 

volume within which BH or NS is to be detected should be proportional M^^^ . Then the ratio 
of detection rates of BH and NS can be written as 

Dbh ^ Rbh /MBj^y^/' 
Dns Rns \ Mns J 

Let us make a simple estimate. Take a GW detector with hj-ms = 10^^"*^ and S/N = 1 (as 
for the initial LIGO or VIRGO interferometer). Let the mass of NS and BH be 1.4M0 (the 
typical value well justified experimentally) and 8.5M0 (the mean mass of the BH candidates, 
see Table Q). Let us also assume that any star with the initial mass M{NS) > lOM© yields 
NS (the typical value confirmed theoretically), and the threshold mass for BH formation is the 
maximal from the estimates given above, M{BH) > Mcr = SOM© . Hence using Salpeter mass 
function we find 

Rbh N{M > SOM©) /SOM© 



Rns N{M>WMq) \WMq 

and then 

Dbh /8OM0\-1-35 / 8.5Mq \ ^^/^ 



-1.35 

~ 0.06 



_ 0.06 X 90 

Dns \WMqJ \1A0Mq; 
This estimate is, of course, very rough, but the precise calculation gives essentially the same 
result. Fig. § shows the calculated absolute registration rate of mergings D by the detector 
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Figure 5: A schematic dependence of the detection rate ratio on the detector sensitivity (the upper scale) and the 
limiting distance from which binary NS (1.4 + 1.4M0, the middle scale) and binary BH (8.5 + 8.5M0, the bottom 
scale) mergings can be detected. See the text for more detail. 



with hrms = 10~^^ and S/N = 1, as a function of kBH- The calculations have been done for 
Lyne-Lorimer kick velocity distribution with the parameter wo = 400 /. The vertical dispersion 
is due to the dependence of D on Mcr- It is seen that binary NS mergings occur with a rate of 
about 1 event per year. In the same time, the total detection rate (except for ksH ~ 0.1) can 
exceed that of binary NS mergings by 2-3 orders. The filled region shows the region of the "most 
probable parameters", for which the calculations meet both the upper limit to the number of 
BH with radiopulsars (less than 1 per 700 single pulsars) and the number of BH candidates of 
Cyg X-1 type (from 1 to 10 in the Galaxy). In this region, from 10 to 100 mergings per year 
must be detected, with the majority of them being binary BH mergings. 



4-1 Discussion 

Nevertheless, one may wonder whether binary BH (which are less numerous than binary NS) 
should always be detected more frequently than binary NS. It seems worth looking more closely 
on the stellar mass distribution around a GW detector on Earth (Fig. |5|). In this Figure, the 
detection rate ratio Dbh / Dj^s is plotted schematically against the detector sensitivity level (the 
upper scale) , which can be expressed through the maximum distance from which binary NS /BH 
mergings with M^s = 1.4M0 and Mbh = S.SMq can be detected (two bottom scales). Four 
segments may be distinguished in this plot (from left to right): first, when we register objects 
inside some part of the Galaxy, second, when all objects within the Galaxy are detected but 
no extragalactic objects can be detected, third, mainly extragalactic events are detected from 
distances more close than those at which the initial star formation occurs, and forth, where 
we detect all events in the Universe. In different segments different detection ratios will be 
obtained. In the first segment, the detection ratio depends on the galactic structure: if NS and 
BH populate the same spherical halo, this ratio (roman HI on the vertical axis; the solid line) is 
Rbh /Rns{Mbh /Mns)^^^^ ; if NS and BH populate the galactic disk, this ratio becomes (roman 
II, the bottom dashed line) Rbh / Rns{Mbh / MnsY^'^ ; if NS fill more extended halo than BH, 
i.e. the halo radius r^s > tbh (roman IV, the upper dashed line), then the detection ratio is 
Rbh /Rns{Mbh /MNsy^^^{rBs/fNs)^ ■ III the second and fourth segments the detection ratio 
will be minimal and simply equal to (roman I) Rbh / Rns ■ In the third segment, the detection 
ratio of type III is realized. At the end of this segment evolutionary effects can affect the 




Figure 6: (Left panel) log A^-log S curves calculated for different spectral power-law indices attributable to gamma- 
and X-ray emission in a GRB. Values of the cosmological model parameters are shown in the Figure. The bar 
indicates the accessible range of the total GRB rates in the entire Universe varying the parameters as discussed 

in the text. 

Figure 7: (Riffikt panel) The 2-nd BATSE catalog (solid points) is fitted by the cosmological GRB model (from 
Lipunov et a/.Ej). Note that the total GRB rate in the Universe is ~10'* per year, ~3 orders of magnitude smaller 

than the total binary NS merging rate. 



detection ratio (the dashed line). 

Note that the present sensitivity of bar detectors {hrms ~ lO^"*^^) falls within the second 
segment (an "unhappy" situation because no mass-ratio enhancement for BH detection occurs), 
whereas the initial laser interferometers {hrms ~ 10~^^) are "luckily" in the third segment with 
the enhanced type III detection ratio. With the advanced LIGO sensitivity, it is possible to 
detect evolutionary effects and even reach the very edge of star formation. 

5 Gamma-Ray Bursts 

Using the obtained above dependence of compact binary merging rates in the elliptical galaxy 
on time (Fig. ||) and assuming the cosmological origin of GRB as products of binary NS/BH 
coalescences, we can compute the theoretical logA^-logS curve. To do this, we need to specify 
the cosmological parameters, the moment of the star formation beginning, and the spectral 
power-law index of a typical GRB (see Lipunov et a/. for more detail). Taking the density 
of baryons in stars 17* = 0.0046 (in terms of critical density to close the Universe) E3 and 
varying other parameters within limits permitted by the present theory and observations (A- 
term: < Q\ < 0.7; the fraction of elliptical galaxies: 0.15 < e < 0.9; the star formation starting 
redshift: 2.5 < < 10), the total compact binary merging rate in the Universe (a constant the 
logA^-log/S curve goes at small fluxes) is found to vary not too much from a few- 10^ yr~^ to a 
few-10^ yr~^ (see Fig. |). 

This curve is consistent with observational data obtained by BATSE (Fig. |^ from0). How- 
ever, the total GRB rate in the Universe will be of order ~10^ yr^^- The difference in about 3 
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Figure 8: The observed log A'^-Iog S distribution for the Einstein IPC soft X-ray flashes (from Gotthelf et al.U; the 
broken line) and theoretical log A^'-log S curves (with the same parameters as in Fig. n) for four different values 
of the mean binary NS merging rate in a spiral 10^^ galaxy. 



orders can be explained either by assuming the only one me 
the gamma-ray emission collimation into a ~7° soHd angle t 



ring of ~1000 to yield GRB, or by 



6 Einstein Soft X-ray Bursts 

In 1996, after Einstein IPC archive data reprocessing, weak soft X-ray flashes were reported!, 
which are distributed isotropically on the sky and unidentified with known astronomical objects, 
like GRB. The interpolation of the Einstein IPC field of view on the total celestial sphere yields 
the rate of these flashes ~ 2 • 10^ yr~^. If the existence of these flashes is confirmed (for example, 
using independent ROSAT and ASCA data) and they are generated by binary NS coalescence, 
they will be a serious indication favouring the high, '^1/10000 yr^^, NS merging rate in the 
Galaxy (see Fig. H). 

7 Conclusion 

To conclude, we expect the GW events due to compact binary coalescences to be registered 
already by the initial laser GW- interferometers of LIGO/VIRGO type at a rate substantially 
higher than has previously been thought. We also think that BH will thus be detected simultane- 
ously with GW. Otherwise, we will have very stringent constraints on BH formation parameters, 
which is, of course, much less interesting. 
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